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We report on recent studies of QCD performed using reanalysed e + e annihilation data 
recorded at centre of mass energies 14 < y/s < 44 GeV by the JADE experiment oper- 
ated 1979 to 1986 at the PETRA e + e" collider at DESY, Hamburg, Germany. The data for 
some event shape observables are compared to modern Monte Carlo event generators using 
parameter values obtained from LEP 1 data. The distribution of £ = ln(l/a;) with x = 2p/^/s 
and p the momentum of a charged hadron is measured and compared to QCD predictions. 
Fits of 0(as)+NLLA (resummed) QCD predictions combined with power corrections to event 
shape data including for the first time the yfs = 14 and 22 GeV data samples are discussed. 



1 Introduction 

QCD as the theory to describe many features of hadron production in e + e~ annihilation such 
as the total cross section, jet formation, multiplicity or momentum spectra has been established 
mainly with data from the experiments at the PETRA e + e~ collider, see e.g. ^I^IIHIIUIHI- 
However, the analyses of that time were limited by the theoretical understanding of QCD as 
well as the models needed to correct for hadronisation effects. 

In particular during the LEP era from 1989 to 2000 enormous progress in improving QCD 
calculations by adding higher order radiative corrections and creating generic numerical integra- 
tion programs jHUZj or resummed logarithmic contributions (e.g. 8J was made. Also, detailed 
predictions directly relating perturbative QCD (pQCD) predictions with hadron distributions 
for e.g. multiplicity or momentum spectra became available, see e.g. |9*UTU]. 



The modelling of hadronisation effects had benefited from the use of Monte Carlo event gen- 
erators such as P YTHIA / JETSET [H], HERWIG [E] or ARIADNE H3j containing consistent 
descriptions of parton and hadron production tuned to precise data from the LEP experiments, 
see e.g. |14|I15| . In another approach analytical models of hadronisation applicable to some event 
shape or jet production observables had been developed allowing comparisons of theory with 
data without Monte Carlo hadronisation models, see e.g. |16j . 

PQCD effects grow with decreasing energy scale Q like 1/lnQ and in particular non- 
perturbative effects in many event shape or jet production observables grow with 1/Q. Therefore 
the data from the JADE experiment can still be used for important and unique tests of modern 
QCD predictions. From another point of view one can argue that only now our understanding 
of QCD and hadronisation processes is sufficiently advanced to allow consistent interpretation 
of the PETRA data. 

The JADE experiment is described in The data used in the studies presented here 
correspond to samples of about 1500 events at = 14, 22 and 38 GeV, about 4000 events at 
44 GeV and about 35000 events at -y/s = 35 GeV. The JADE software consisting of a detailed 
detector simulation program and the reconstruction programs has been installed on an IBM 
RS6000 AIX system using the IBM Fortran compiler |17j . It is possible to generate final states 
from e + e~ annihilation using modern generator programs like PYTHIA/ JETSET, HERWIG or 
ARIADNE, to pass these through the JADE detector simulation and to finally reconstruct the 
output of the detector simulation in the same way as the data. 

2 Data vs modern MCs 

Before the modern Monte Carlo generators can be used with confidence in the data analysis 
they must be compared with data to establish their performance. The data for the event shape 
observable thrust T and the distributions of y23 a shown in figure ^ are corrected for acceptance, 
detector resolution and initial state radiation (ISR). In addition the contribution from e + e~ — > 
bb events has been subtracted using simulated events from PYTHIA. The e + e~ — ► bb events, 
approximately 9% of the JADE samples, are considered as background since the high masses 
of the b-hadrons are known to distort e.g. event shape distributions, in particular at lower 
\fs |18U19llT7| . The data are compared with predictions based on u, d, s and c quarks from 
the generator programs PYTHIA 5.7, ARIADNE 4.08, HERWIG 5.9 and COJETS 6.23 as well 
as JETSET 6.3 as used in earlier JADE analyses. The parameter sets for all programs except 
JETSET are taken from the OPAL collaboration, i.e. they are derived from LEP 1 data at 
-^/s ~ M z o, while the JETSET parameters are taken from JADE, see for details. 

We find that the PYTHIA and ARIADNE generators generally describe the JADE data 
well. HERWIG is seen to produce somewhat too many 3-jet like events and correspondingly too 
few 2-jet like events at lower values of y/s while COJETS does not give a reasonable descrip- 
tion of the data. The COJETS model implements a leading-log parton shower without taking 
colour coherence into account combined with independent fragmentation j2Uj. The JETSET 6.3 
predictions also describe the data fairly well. In conclusion we find that modern Monte Carlo 
Generators using parameter sets derived from LEP 1 data can be used for QCD analyses with 
JADE data. 

3 Scaled Momentum Distribution 

The shape of the distribution of £ = ln(l/x) with x being the momentum of charged particles 
normalised to the beam energy is predicted in pQCD to be approximately a skewed Gaussian, see 

"The value of the jet resolution parameter where the event changes from a 2-jet to a 3-jet configuration with 
the Durham jet algorithm. 




Figure 1: Comparison of JADE data for 1-T (left) and j/23 (right) corrected for experimental and ISR effects and 

bb-background with Monte Carlo predictions |17) . 



e.g. |9*HTU] , The position of the peak £° is predicted as a function of The pQCD predictions 
in the modified leading log approximation (MLLA) including subleading effects |21lll0ll2"2] may 
be directly compared with the hadron-level data up to a normalisation factor, a property known 
as local parton hadron duality (LPHD). 

The £ distributions have been measured using JADE data at 22, 35 and 44 GeV [22]. In 
this analysis no correction for the presence of e + e~ — ► bb events has been made consistent with 
analyses by other experiments. Figure (left) shows the result for y/s = 35 GeV. Superimposed 
is a fit of the pQCD prediction to determine the position of the peak £°. The figures for the 
other y/s points at 22 and 44 GeV look similar but have larger statistical errors due to the 
smaller data samples. 

The fit is seen to describe the data well around the peak region; the prediction is not expected 
to be in good agreement with the data outside of the fit region. The results for the peak positions 
and effective QCD scales are given in table ^ 



Table 1: Results for peak position £° and A c ff from pQCD fits to JADE data with total uncertainties |22|. 



[GeV] 


£ U 


Acff[MeV] 


22 


2.74 ± 0.09 


136 ± 28 


35 


3.06 ±0.05 


142 ± 25 


44 


3.19 ±0.06 


110 ±38 




Figure 2: The figure on the left shows the distribution of £ measured at y/s = 35 GeV corrected for experimental 
effects. The figure in the middle shows the JADE results for £° together with results derived from OPAL data 
and compared with a pQCD fit for the */s dependence of £°. The figure on the right shows the same results for 
f° except at i/s ~ M z o where results from flavour separated data from OPAL are given; the fits are explained in 

the text 



The dependence on yfs of £° is shown in figure |2] (middle) . The three JADE data points are 
compared with results derived using the same fits from OPAL data and fitted with the pQCD 
prediction. The fit describes the data reasonably well, however, the result for A e ff differs from 
the results shown in tabled by two to three standard deviations. 

In order to study effects of heavy flavours in the data samples a fit with an effective QCD 
coupling A e ff for each flavour (uds, c or b) is made. The three A e ff are related by £ Cj b = Cuds + 
0.5 ln(A c b/A uc js) and flavour separated data from OPAL are used at yfs ~ M z o. At the other 
yfs the known flavour fractions in e + e~ — ► qq are employed. The results A u d s = 184 ± 33 MeV, 
A c = 239 ±90 MeV and Ab = 247 ±28 MeV indicate the presence of quark mass effects of about 
20 to 30%; the fits are shown in |2] (right). 

4 Extended Power Correction Fits 

In this section we return to the study of event shape and jet production observables, describing 
results of ^7]. We concentrate on fits of pQCD predictions in O(o^) matched with resummed 
NLLA calculations as described e.g. in |19U16ll23j . In these fits the hadronisation corrections are 
implemented using the analytic DMW model • The strong coupling ag is assumed to remain 
finite at very low scales below and above the Landau pole such that it can be integrated: ao(fii) = 
Jq 1 as(k)dk with m = 2 GeV. The value of ao cannot be predicted and is a free parameter 
of the model. The main prediction of the model for differential event shape distributions of some 
selected observables is that a shifted pQCD prediction describes the corresponding hadron-level 
data: da/dy = (da/dy)pT(y — D y P), where D y is an observable specific factor or function and 
P is an universal factor: P ~ fj,i/Q(ao — as). Thus P contains all non-perturbative dependence 
on the scale Q = i/i and the parameter ao- 

We want to discuss fits of such predictions to the observables Bw and 2/23 which have special 
properties with standard fits. With Bw global fits to data from ^/s = 189 to 14 GeV are found 
to describe the data less well compared to other observables, in particular at low -y/i [TT|IT9]. 
Attempts were made to add extra terms to the power correction part of the prediction and a 
term of the form An ln(Q) /Q was found to provide the best improvement, albeit with a negative 
coefficient An = — 0.60 ± 0.06(fit). The x 2 /d.o.f. value of the fit improved from 134/154 for the 
standard fit to 77/153. The results are shown in figure (left), displaying for clarity only data 
at y/s < Mzo. The dotted lines are the standard pQCD combined with power correction fits 
which are seen to deviate from the data in particular at y/s = 14 and 22 GeV. The fit of the 



same prediction with a An ln(Q)/Q term added shown by the solid lines improves the agreement 
significantly. There is thus evidence for the presence of additional terms in the By/ prediction 
probably behaving like ln(Q)/Q. 



l/o do/dB l/o dc/dy 



o 




Figure 3: The figures show data at y/s = 91 to 14 GeV for Bw (left) and j/23 (right), including data from JADE. 
The data are compared with pQCD combined with standard power corrections and modified power corrections 
as indicated on the figures |17|. The solid lines indicate the fit ranges for all fit variations. 

The DMW prediction for (1J23), the mean value of the distribution, is that its power correction 
is proportional to ln(Q)/Q 2 24 . Previous analysis of mean values have supported a strongly 
suppressed power correction ^F3\. Here the new data for j/23 from JADE are used to study the 
?/23 distribution [23^3 ■ Figure |3] (right) shows fits of pure pQCD to the data as dotted finest 
The fits are seen to describe the data well at large y/a but deviations become visible within the 
fit ranges at low yfs = 14 and 22 GeV. The fit of pQCD with a A20/Q 2 term added is shown 
by the solid and dashed lines. These fits are indistinguishable from the pure pQCD fits at large 
\fs but have improved agreement at low \fs. The x 2 /d.o.f. values changed from 151/107 for 
the pure pQCD fit to 71/106 for the fit with the additional A20/Q 2 term. The parameter is 
found to be ^20 = 2.25 ± 0.18(fit) GeV 2 . The new data at 14 and 22 GeV play an important 
role in discriminating between the two possible predictions. We thus find evidence for a power 
correction probably behaving like 1/Q 2 at low yfs for 2/23 distributions. 

5 Summary 

The continuing analysis of the e + e~ annihilation data of the JADE experiment provides unique 
opportunities for tests of QCD. Modern Monte Carlo generators with parameter sets derived 
from LEP 1 data are found to describe event shape and jet production observable distributions 
reasonably well, with PYTHIA and ARIADNE showing a somewhat better agreement at low 
y/s compared to HERWIG. New measurements of the charged particle momentum fractions £ = 
ln(l/x) have been presented and compared with MLLA pQCD predictions including subleading 
effects. Studies of extended DMW power correction fits using the new JADE data at ^/s = 14 
and 22 GeV have provided evidence for additional power correction like terms in the DMW 
power correction for the observables By/ and j/23 ■ 

b We don't yet use the improved numerically resummed prediction |2(il 
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